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Abstract 

We investigate the classical limit of the Knizhnik-Zamolodchikov-Bernard equa- 
tions, considered as a system of non-stationar Schroodinger equations on singu- 
lar curves, where times are the moduli of curves. It has a form of reduced non- 
autonomous hamiltonian systems which include as particular examples the Schlesinger 
equations, Painleve VI equation and their generalizations. In general case, they are 
defined as hierarchies of isomonodromic deformations (HID) with respect to chang- 
ing the moduli of underling curves. HID are accompanying with the Whitham 
hierarchies. The phase space of HID is the space of flat connections of G bundles 
with some additional data in the marked points. HID can be derived from some 
free field theory by the hamiltonian reduction under the action of the gauge symme- 
tries and subsequent factorization with respect to diffeomorphisms of curve. This 
approach allows to define the Lax equations associated with HID and the linear 
system whose isomonodromic deformations are provided by HID. In addition, it 
leads to description of solutions of HID by the projection method. In some special 
limit HID convert into the Hitchin systems. In particular, for SL(A^, C) bundles 
over elliptic curves with a marked point we obtain in this limit the elliptic Calogero 
A^-body system. 



1 Introduction 

1.1 It is known for a long time that equations provided isomonodromic deformations 
posses structures typical for integrable systems, such as the Lax representation 0. 
This paper concerns the studies of isomonodromic deformations in the spirit of Hitchin 
approach to integrable systems 0. Hitchin discovered a wide class of hamiltonian in- 
tegrable systems on the cotangent bundles to the moduli space of holomorphic vector 
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bundles over Riemann curves. Some important facts about the Hitchin systems became 
clear later. 

First of all, the Hitchin systems are related to the Knizhnik-Zamolodchikov- Bernard 
(KZB) equations Q on the critical level. KZB equations on the critical level take the 
form of second order differential equations, and their solutions are partition functions 
of the Wess-Zumino-Witten (WZW) theory on the corresponding Riemann curves. It 
turns out that these operators coincide with quantum second order Hitchin Hamiltonians 
@) 0' B Uni' 0]- KZB eqs. for correlators of vertex operators need to include in the 



construction curves with marked points, where vertex operators are located. The KZB 
Hamiltonians for correlators define very important classes of quantum integrable equations 
such as the Gaudin equations (for genus zero curves), elliptic Calogero equations (for 
genus one curves), and their generalizations fl^, [l^, It means that their classical 
counterparts are particular cases of Hitchin systems [0, ^. 

Hitchin approach, based on the hamiltonian reduction of some free field hamiltonian 
theory, claims to be universal in description of integrable systems. Essentially, it allows 
to present almost exhaustive information - integrals of motion. Lax pairs, action-angle 
variables, explicit solutions via the projective method. 

The main goal of this paper is to go beyond the critical level in the classical limit of 
KZB eqs. and repeat the Hitchin program as far as possible. For generic values of level 
the KZB equations have the form of non-stationer Schrodinger equations, where the role 
of times is played by the coordinates of tangent vectors to the moduli space of curves 



1^, 10, 15 1 . On the classical level they correspond to the non-autonomous hamiltonian 
systems. We will call them the hierarchies of isomonodromic deformations (HID). In this 
situation the analog of the Hitchin phase space is the moduli space of flat connections 
A over Riemann curves „ of genus g with n marked points. While the flatness is the 
topological property of bundles, the polarization of connections A = {A, A) depends on 
the choice of complex structure on We consider a bundle V over the moduli space 
Aig,n of curves with fiat connections {A, A) as fibers. The fibers is supplemented by 
elements of coadjoint orbits Oa in the marked points Xa- The points of fibers {A, A) are 
analogs of momenta and coordinates, while the base Aig^n serves as a set of "times" . There 
exists a closed degenerated two-form u on V, which is non degenerated on the fibers. The 
connection A play the same role as in the Hitchin construction, while A replaces the Higgs 
field. Essentially, our construction is local - we work over a vicinity of some fixed curve 
Sg^„ in M.g,n- As we already have mentioned, the coordinates of tangent vector to A4g,n 
at „ play role of times, while in the Hitchin times have nothing to do with the moduli 
space. The Hamiltonians of HID are the same quadratic Hitchin Hamiltonians, but now 
they are time depending. 

There is some free parameter k (the level) in our construction. On the critical level 
(k = 0), after rescaling the times, HID convert into the Hitchin systems. As the later, they 
can be derived by the symplectic reduction from the infinite affine space of the connections 
and the Beltrami differentials with respect of gauge action on the connection. In addition, 
to come to the moduli space we need the subsequent factorization under the action of the 
diffeomorphisms of „, which effectively acts on the Beltrami differentials only. Apart 
from the last step, our approach is closed to |TB[, where the KZB systems is derived as 
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a quantization of the very similar symplectic quotient^ Due to this derivation, we find 
immediately the Lax pairs, prove that the equations of motion are consistency conditions 
of the isomonodromic deformations of the auxiliary linear problem, and, therefore, justify 
the notion HID. Moreover, we describe solutions via linear procedures (the projection 
method). 

For genus zero our procedure leads to the Schlesinger equations. This case was dis- 
cussed earlier . We restrict ourselves to simplest cases in which we consider only simple 
poles of connections. Therefore, we don't include in the phase space the Stokes param- 
eters. This phenomenon was investigated in the rational case in detail in [Q. Isomon- 
odromic deformations on genus one curves were considered in and on higher genus 
curves in Here we consider genus one curves with one marked point and obtain for 
SL{2, C) bundles a particular family of the Painleve VI equations. Generalization of this 
case on arbitrary simple groups leads to the multicomponent Painleve VI related to these 
groups. If we introduce a few marked points we come to the elliptic generalization of the 
Schlesinger equations. In fact, the concrete systems we consider here are the deformation 
from the critical level of those Hitchin systems that described in 0]. 

1.2. Painleve VI and Calogero equations The very instructive (but not generic) 
example of our systems is the Painleve VI equation. It depends on four free parameters 
{PV 10,13,^,8) and has the form 

d^X _1 1 1 1 , ,c/X 2 .1 1 1 JX 

'W ~ 2^x ^ X - 1 ^ x-r^~dF^ ~H^t-i ^ x-r~di^ 

X{X-l){X-t). t-l ,t{t-l), 

PVIa,i3,j,s has a lot of different applications (see, for example |22[). It is a hamiltonian 



systems |23]. We will write the symplectic form and the Hamiltonian below in another 



variables. This equation was derived firstly by Gambler p4|, as an equations which has 



not solutions with moveable singularities. Among distinguish features of this equation 
we are interesting in its relation to the isomonodromic deformations of linear differential 
equations. This approach was investigated by Fuchs p5|. 



There exsits elliptic form of PVIa^p^^^s derived by Painleve itself p6|. It was inves- 



tigated recently by Manin |^ in connection with Frobenius varieties. It sheds light on 
connection of the Painleve equations with the Hichin systems, and, thereby, with the KZB 
equations. We present shortly this approach. 

Let p{u\t) be the Weierstrass function on the elliptic curve = C/(Z + Zr), and 

ei = P(| k), (To, . . . , T3) = (0, 1, r, 1 + r). 
Consider instead of (X, t) in (|1 . 1| ) the new variables 

(^,r)^(X= ^(^'"^-^\ t = gi^). (1.2) 

62 - ei 62 - 61 



^Quantization of isomonodromic deformations on rational and elliptic curves and their relations to 
KZB was considered in H, ^ 
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Then PVIa,i3,^,5 takes the form 



— = -duU{u\T), 
U{u\r) = -^Y.^a,p{u + ^\T), (1.3) 
(ao, . . . , as) = (a, -P, 1,^-^)- 



The hamiltonian form of ( |1.3|) is defined by the standard symplectic form 

uo = 6v6u, (1.4) 

and the Hamiltonian 

H = — + U{u\r). (1.5) 



The equation of motion ( |1.3| ) can be derived from the action 5* 

6S = v6u- H6t. (1.6) 

We can consider (0), (|1.5| ) as non-autonomous hamiltonian system with 
the time-dependent potential. 

To find symmetries we consider two-form 

u = uq — 5H5t = 6v6u — 6H6t. (1.7) 

The semidirect product of Z+Zr and the modular group acting on the dynamical variables 
(f , u, t) are the symmetry of ( |1.7|) . We consider them in detail in Sect. 7. 
Let us introduce the new parameter k and instead of (|1.6| ) consider 



u! = luq 5H5t. (l.t 



Then, (|5.1| ) takes the form 



'^'^ = -duUiu\T). (1.9) 

It corresponds to the overall rescaling of constants aj — >■ Put r = tq + nt^ and 
consider the system in the limit k —>■ 0. We come to the equation 

-d^U{u\To), (1.10) 



{dt^y 

It is just the rank one elliptic Calogero-Inozemtsev equation P5| , 3D|, which we denote 
CIa,i3,'y,s- Thus, wc have in this limit 

PVIa,l3r(,S '^—^ C 10,13,^,5 ■ (1-11) 
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So far we don't now how to manage with the general forms of both types of equations 
for arbitrary values of constants. Here we consider only the one-parametric family in (|1.9|) 

4 ' 4 ' 4 '2 4 

J 4'^ 4'^ 4'' 4' 2 4 
The potential (|1.9|) takes the form 

^(^1^) = 7^^V(2«|r), (1.12) 



We will prove that ( |1.4|) (|1.5| ) with the potential ( |1.12|) describe the dynamic of fiat 



connections of SL(2, C) bundles over elliptic curves Tr with one marked point Si i. In 
fact, u lies on the Jacobian of T,-, (f , u) are related to the fiat connections, and r defines a 
point in Roughly speaking, the triple (f , m, r) are the coordinates in the total space 

of the bundle V over A^i,i, uj ( |1.8|) is the two- form on P, and cuq is the symplectic form on 
its fibers. This finite-dimensional Hamiltonian system is derived as a quotient of infinite- 
dimensional phase space of {A, A) connections and the Beltrami differentials as times 
under the action of gauge transforms and diffeomorphisms of Si i. This approach directly 
leads to the Lax linear system and allows to define solutions of the Cauchy problem via 
the projection procedure. Simultaneously, we describe the auxiliary linear problem whose 
isomonodromic deformations are governed by this particular family of the Painleve VI. 
The discrete symmetries of (|1.6| ) are nothing else as the remnant gauge symmetries. On 
the critical level it is just two-body elhptic Calogero system. The corresponding quantum 
system is identified with the KZB equation for the one-vertex correlator. In the similar 

way PVI^ is the classical limit of the KZB for /t 7^ 0. This example will be 

4 ' 4 ' 4 ' 2 4 
analyzed in detailed in Sect. 7. 



Solutions of PVI^ _^ ^ in particular case (z/^ = ^) were found by Hitchin [RI 
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in connection with his investigations of anti-self-dual Einstein metrics. They are written 
in terms of theta functions. But for generic values of u solutions are genuine Painleve 



transcendents |^ . 



1.3. Whitham equations The previous example is not entirely exhaustive. A special 
phenomena occurs if dimAig n > 1. We have as many Hamiltonians as dimAig n - each 
Hamiltonian Hg is attached to the tangent vector ts to J^g,n at some fixed point 
There are consistency conditions of the equations of motion for the non-autonomous multi- 
time hamiltonian systems. They take a form of the classical zero-curvature conditions for 
the connections 

dt,+H,, (1.13) 

where the commutator is replaced by the Poisson brackets. The later is the inverse tensor 
to the symplectic form, defined on the fibers. The flatness conditions is the so-called 
the Whitham hierarchy, (WH) defined in the same terms by Krichever though the 



starting point of his construction is different and based on the averaging procedure. Let 
S be the action for HID {SS = 6^^uj as in( |1.6|) ). The r-function of WH is 

loe; T = S. 
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It allows to find the Hamiltonians. For HID in the rational case the r-functions were 
investigated in 0. For the Painleve I- VI they were considered in p4| . 

The quantum analog of the Whitham equations was exploited in ||15|, ^ to construct 



the KZB connections. The quantum version of ( p..l3| ) is an operator, acting in the space 



V of sections of holomorphic line bundle over moduli of fiat connections. The quantum 
WH is the flatness of the bundle of projective spaces PV^ over the Teichmiiller space 7^ „. 

1.4. Outline We start in sect. 2 with a general setup about fiat bundles over singular 
curves. In sect. 3 we present the basic facts about the abstract nonautonomous hamil- 
tonian equations. They defined by a degenerate closed two-form, on the extended phase 
space, which is a bundle over the space of times. The Whitham equations occur on this 
stage. We remind the symplectic reductions technique, which is applicable in degenerate 
case as well. Sect. 4 contains our main result - the derivation of the HID, corresponding 
to the fiat bundles over Riemann curves. In sect. 5 we discuss two limits - the level zero 
limit of HID to the Hitchin systems and the classical limit of general KZB equations to 
HID. Then we analyze in detail two feasible examples - fiat bundles over rational curves, 
leading to the Schlesinger equations (sect. 6), and flat bundles over elliptic curves, re- 
sponsible for the Painleve VI type equations, and the elliptic Schlesinger equations (sect. 
7). In particular, we apply the projection method to construct perturbatively solutions of 
the Schlesinger equations. In the Appendix we summarize some results about the elliptic 
functions, which are used in sect. 7. 
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2 Flat bundles over singular curves 

We will describe the general setup more or less naively. We don't consider in this section 
the mechanism of symplectic reduction in detail and postpone it on Sect. 4 We will consider 
three cases: 1) smooth proper (compact) algebraic curves; 2) smooth proper algebraic 
curves with punctures; 3) proper algebraic curves with nodal singularities (double points). 

Let G be a semisimple group and V be its exact representation; f.e. G = SL(A^, C) 
and V is a A^-dimensional vector space with a volume form. 

2.1. Smooth curves. 

Let S" be a smooth oriented compact surface of genus g. Let us consider the moduli 
space FBuns^G of flat V-bundles on S. This space can be considered as the quotient of 
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the space FConn of flat C°° connections on trivial ^-bundle by action of the gauge group 
Q of G- valued C°°-functions on S, or as a result of hamiltonian reduction of space Conn 
of all connections by action of the gauge group. The symplectic form on the space Conn 
is the form 

u= [ <6A,6A >, (2.1) 
Js 

where < •, • > denotes the Killing form, 6A is a Lie(G)-valued one-form on S. So, 
< SA, 6 A > is a two-form, and the integral is well defined. The dual space to the 
gauge algebra Lie is the space of Lie((j)-valued two-forms on S, and the momentum map 
corresponds to any connection A its curvature F4 = dA + \ [A, A\ . Hence, the preimage 
of under momentum map is the space of fiat connections. 

A complex structure S on 5* is a differential operator d : fi^oo fi^oo; its kernel is the 
space of holomorphic functions. For any connection we can consider 9-part of connection; 
it defines a holomorphic bundle. A section is holomorphic, if it belongs to the kernel of 
this operator. Let us recall that: 

1. Connection on the holomorphic bundle £ is the operator of type d + A. 

2. Connection is holomorphic if this operator commutes with operator of complex struc- 
ture on the bundle. 

3. Connection is fiat if its square vanishes; for curves this condition is empty. 

4. The space of all holomorphic connections on £ is an infinite dimensional affine space 
and the gauge group Q acts on it by affine transformations. 

5. The quotient of stable holomorphic connections with respect to the gauge transforma- 
tions is the moduli space Bunj^^c of holomorphic vector bundles on S. 

Let 6 be a point in BunY,,G and Vb is the corresponding V^-bundle on S. Denote by 
adh the bundle of endomorphisms of Vb as bundle of Lie algebras and by Flatb the space 
of flat holomorphic connections on holomorphic bundle Vb- 

Consider the map from FBuns^c onto the moduli space BunY,,G- The fiber of the 
projection FBuns^c ~^ Bunj^^c over a point b is naturally isomorphic to Flatb- These 
fibers are Lagrangian with respect to u - any complex structure on S defines a 



polarization on FBuns,G- 

The tangent space to BunY.,G at the point h is canonically isomorphic to the first 
cohomology group H^{ll,adb)- The space Flatb of holomorphic connections on Vb is an 
affine space over the vector space adb ® of holomorphic ac?f,-valued one-forms, 

since a difference between any connections is a ac?f,-valued differential form. The vector 
spaces adb) and adb ® are dual. 

The Dolbeault representation of cohomology classes corresponds to the description 
based on the hamiltonian reduction. Indeed, let us decompose A into (1,0) and (0,1) 
parts: A = A + A. From the zero curvature condition we get that (1, 0)-part 5 A must be 
A-holomorphic: Sa^A = [d + A)5A = 0, and 6 A is defined up to the infinitesimal gauge 
transformations 6 A — > 6 A + dAh. 

2.2. Curves with punctures 

For any two isomorphic representations V and V of G denote by Isom{V, V) the 
space of G-isomorphisms between them; this space is a principal homogeneous space over 
G- For a curve S„ with n marked points xj we replace the moduli space Bun-^^c by moduli 
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space BunY,,yi,G-i (x = (a;i, . . . of holomorphic ^-bundles V with the triviahzations 
9] '■ — > y of fibers at the marked points. We have natural "forgetting" projection 
TT : Bun-s^^^G Buns^G- The fiber of this projection is the product H Isom{Vxj-i V) of the 
spaces of isomorphisms between fibers of bundle at marked points and V . The projection 
TT can be treated as a reduction by action of n copies of the group G which acts transitively 
on the fibers 

BunY„G = BunY:,^^G/ H 

The bundle of endomorphisms of this data is the bundle a(i(,(— x) of endomorphisms, 
vanishing at the marked points, since non vanishing endomorphisms change triviahzations. 
Hence, the tangent space to -Biins,x,G is isomorphic to H^{T,,adb{—x)). The dual space 
to this space is adi){y:) <S)fl^). Consequently, in order to get the symplcctic variety 

we must replace the affine space Flath of fiat holomorphic connection by the affine space 
Flath{log) of fiat connections with logarithmic singularities. As a result, we get the 
moduli space FBunY,,yi,G of triples: (holomorphic bundle, triviahzations at marked points, 
holomorphic connection with logarithmic singularities). Again, we have the Lagrangian 
projection Fi?Mns,x,G Bun-E^:x:,G- The fiber of this projection is an affine space over 
H°{^,adb{x)0n^).' 

The product of n copies of G acts on FBunx,:>i,G by changing trivialization of fibers, 
this action is hamiltonian. We can consider the hamiltonian reduction with respect to this 
action. As a result, for any collection of G-orbits of adjoint action we get the moduli space 
FBun-s^:^ Q.{o } of pairs (holomorphic bundle, holomorphic connection with logarithmic 
singularities with residues in orbits Oj). 

According to the Dolbeault theorem, the tangent space H^{J2,adb{—^)) can be re- 
alized as the space of (0, 1) ad^-valued forms holomorphically vanishing at the marked 
points modulo ^^-coboundaries of holomorphically vanishing at marked points ad;,- valued 
functions. A function (or one-form) is referred to be holomorphically vanishing at the 
point X, if it has asymptotics {z — x)0{l) at this point. 

The second description: the tangent space if^(E, a(i6(— x)) is isomorphic to the space 
of all (0, 1) adfe- valued forms modulo 9yi-coboundaries of vanishing at marked points adh- 
valued functions. Indeed, we have a natural embedding of the space of holomorphically 
vanishing forms to the space of all forms and the space of holomorphically vanishing 
functions to the space of vanishing functions. It defines the map from the space from 
first description to space from second description. We can assume that A vanishes at small 
vicinities of marked points. Then equals d at these vicinities. Let h he a function, 
vanishing at the marked points which is not holomorphically vanishing, 

oo ^, 

/i = ^ ai{z — Xj) + {z — Xj)o{l). 

i=l 

Then 

dh ^ (^iai{z — Xj) + {z — Xj)o{l))dz 

i=l 

is not holomorphically vanishes. So, the map is injective. At the other hand, for any 
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form u with asymptotics 

oo 

v ~ ai{z — Xj) + {z — Xj)o{l))dz 

i=0 

at the points Xj we can choose some function (p such that u — has asymptotics 

{z — Xj)o{l), since for any collection of asymptotics {J2^o j^^^ii^ ~ ^j) + (-^ ~ ^j)o(l)) 
a function with such asymptotic exists. So, the map is surjective. 

The third equivalent description of this space is defined as the cokernel of the map: 

Q^l,{adb) flc^^{adb) + ^adb\x.; {dAh,h{xi)). 

Indeed, for any element {SA, {ai}) in the cokernel we can choose its representative with 
vanishing second part {flj}. This choice is imique up to for vanishing at the marked 
points functions h. This description is adapted to the forgetting map tt. The local part 
{tti} corresponds to the tangent space of fiber and the global part 5 A corresponds to the 
tangent space of base. 

By "integration" of the action of the gauge Lie algebra adb to the action of the gauge 
Lie group Q we get two description of the moduli space. 

A. The moduli space i?Mns,x,G is the quotient of the space of 5- connections A on the 
trivial F-bundle by action of the reduced gauge group of G- valued functions, whose values 
at the marked points are equal to the neutral element 1 of the group G. The isomorphisms 
of fibers with V are identical maps. 

B. The moduli space BunY.^yi,G is the quotient of the space of pairs (a (9-connection A 
on the trivial V^-bundle, a collection gj of elements of G) by action of the gauge group Q. 
The isomorphisms of fibers with V are Qi. The isomorphism of these descriptions can be 
proved by the following consideration. For any collection [A, gj) we can choose the gauge 
equivalent collection with gj — 1, and such collections are equivalent up to action of the 
reduced gauge group. 

In what follows we assume the second description of BunxT.,x,G 

Let us consider the space of data ((9-connections A on the trivial V"-bundle, a collection 
gj of elements of G, 9-connection A with logarithmic singularities at marked points). This 
space is symplectic with the form FBunx,{xj},G equals to Hamiltonian reduction of this 
space by action of gauge group Q. 

It is worthwhile to note that (in contrast with Case 1) this construction is essen- 
tially based on the complex structure on the surface, since A is nonsingular and A has 
singularities at the marked points. We can consider the moduli space of fiat bundles 
over noncompact surface S \\Jxj — G- representations of fundamental group oi S \ \Jxj. 
According to Deligne, for any complex structure on S any stable representation can be 
realized by connections with logarithmic singularities on holomorphic bundle. Hence, any 
connection A is gauge equivalent to a connection with regular 5-part, but corresponding 
gauge transformation has any singularities at marked points a priori. This fact makes 
such approach very complicated. 

2.3. Curves with double points 
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A curve with a double point (the nodal singularity) can be treated as a "limit" of 
nonsingular curves under pinching of some circle. If this circle is homological to zero, then 
the resulting curve is the union of two intersecting smooth curves; and sum of genera of 
these curves is equal to the genus of nonsingular curves. If this circle is homologically 
nontrivial, then the singular curve is a smooth curve with "glued" two different points. 
The genus of this curve is less by one that genus of initial curves. The normalization of 
singular curve ( "disglueing" of singularity) is a smooth curve (not connected, in general) 
with marked points. 

Let us fix some notations. Denote by Aig the moduli space of smooth curves of genus 
g, and denote by Aig^n the moduli space of smooth curves of genus g with n different 
marked points {M.g = Aigfi) ■ Let M.g,n be the Deligne-Mumford compactification of 
M.g,n- Then the compactification divisor Dec = M.g \ M.g is the union of components 
D^'^'^, 91+92 — 9 and D^^. These components are covered by A4gi,i xM.g2,i and M.g-1,2 
corresp ondingly. 

Consider the moduli space BuUg j^ of pairs (smooth curve of genus g, G-bimdle on it). 
Evidently, this space is fibered over J\4g with Bunx,G as fibers. Unfortunately, we don't 
now the canonical compactification BuUg^x, which is fibered over A4g. Let us assume 
that such compactification does exist. Then the open part of fiber over singular curve 
must be described as moduli space of the following data (G-bundle over normalized curve, 
isomorphisms of fibers over "glued" points). 

Denote by S° a singular curve with nodal points 1/2 • • • |/n, and denote by S its 
normalization, Xa and are preimages of under normalization map. Then moduli 
space Bun-so Q, corresponding to E° is the moduli space of (holomorphic ^-bundle on 
E, isomorphisms between fibers of this bundle over Xa and Xn+a)- This space is quo- 
tient of BunY,,yi,G by action of n-copies of group G, j-th G acts on /som(Vj;„, V^) and 
Isom{yxn+ai ^) from the right: 

{ha\l<a<n '■ {9a, 9a+n}l<a<n {9aha, 9a+nha}l<a<n- 

Corresponding symplectic variety FBun-so Q is the result of symplectic reduction of 

FBun-E^x,G by the described above action with zero level of momentum map. This space 
is the moduli space of data (holomorphic K-bundle on X, isomorphisms between fibers of 
this bundle over Xa and Xn+a, connection with logarithmic singularities such that residue 
at points Xa and Xn+a are opposite). In the last part of data we use isomorphisms of fibers 
from the second part of data. 

3 Hamiltonian formalism 

Wc consider nonautonomous hamiltonian systems. For this type of systems it is the cus- 
tom in the classical mechanics to deal with a degenerate symplectic form on the extended 
phase space, which includes beside the usual coordinates and momenta, the space of times 
and corresponding to them Hamiltonians as the conjugate variables. Then the hamilto- 
nian equations of motion are defined as variations of dynamical variables along the null 
leafs of this symplectic form. 

3.1. Equations of motion 
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Let TZhe a (infinite-dimensional) phase space endowed with the non-degenerate sym- 
plectic structure. For simphcity, we take it in the canonical form 

uo = (5v,5u), 

V = {Vi, ...,Vi...), U = . . . , Mi . . .). 

Consider the space of "times" Af = {t = (ti, . . . , t^, . . .)} and corresponding dual to 
times Hamiltonians {H^, . . . , if", . . .) on 7^ depending on times as well. We consider the 
extended phase space which is the bundle V over A/" with fibers TZ. Introduce a symplectic 
form on V 

uj = ujQ-Y, dH^dta = (5v, 6u) - J2 SH^6ta. (3.1) 

a a 

This form is closed but degenerated on V. The vector fields 



lie in the kernel of u iff 



Then the vector field V° G Ker uo takes the form 

It can be checked immediately that they commute. For any functions fiV, U,t) on V its 
evolution is defined as 

"^^^^'"'^^ = VV(v, u, t), e Ker u. (3.3) 

dta 

It follows from ( p. 21 ) that the Hamiltonians subject to the classical zero curvature condi- 
tions {the generalized Whithem hierarchy) 

- +{H\H'^U = 0, a,b = l... (3.4) 
atb ata 

Evidently, for the time independent Hamiltonians ( |3.3| ),( p^ give the standard approach. 
Define the action S of the system as S~^uj = 6S 

6S=iY,6u)-J2HJta (3.5) 

a 

and the r-function 

S\ogT = 5S. (3.6) 
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The equations of motion can be written down in the Hamihon-Jacobi form 

dS 6S 

- = -i7.(-,u,t). 

3.2. Moment map 

The symplectic reduction can be apphed to the nonautonomous system with symme- 
tries in the standard way. We shortly repeat this approach. Let ^ be a symmetry group 
of the system. It means that for any e GLie(^) there exists the vector field acting on V 
such that the Lie derivative annihilates u 

C,uj = {5j^ + jj)uj = 0. 

The vector field is called a hamiltonian vector field with respect to uj. Since u is closed 
we can write locally j^u = 6F^. If j^uj is exact then is the strictly hamiltonian vector 
field. The function = F^{\, u, t) is a linear function on Lie(^) 

=< e, J^(v,u,t) > 

and thereby defines the the moment map 

J:V^ Lie*{g). 

Assume that we put the moment constraint 

J(v,u,t)=0. (3.7) 

The quotient 

is a symplectic space with the symplectic form uj^^'^ which is the reduction of the original 
form u (13.11). It is defined by the two step procedure. 

i) Fixing the gauge. In other words, there should be defined a surface in the phase space 
V, which is transversal to the orbits of the gauge group Q. 

ii) Solving the moment constraint equations ( p.7|) for the dynamical variables, restricted 
to the gauge fixing surface. 

The main difference with the autonomous hamiltonian systems is that we demand the 
gauge invariance for the whole degenerate symplectic form (p.l|) only and don't consider 
the form ujq and the Hamiltonians separately. 

4 Symplectic reduction and factorization 

4.1. Space of "times" Af' . 

Let Tjg^n be a Riemann curve of genus g with n marked points . Let us fix the complex 
structure on Eg^„ by defining local coordinates {z, z) in open maps covering T^g^n- Assume 
that the marked points (xi, . . . ,x„) are in the generic position, i.e. there exists a set of 
their vicinities (Wi, . . . ,Z^n), such that fl = for a 7^ h. 
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The deformations of the basic complex structure are determined by the of Beltrami 
differentials /i, which are smooth (—1, 1) differentials on We identify this set with the 
space of times M' . The Beltrami differentials can be defined in the following way. Consider 
a chiral smooth transformation of Sg,n; which in some local map can be represented as 

w = z — e{z,z), w = z (4.1) 

Up to the conformal factor 1 — de{z, z) the corresponding one-form dw is equal 

daiz z) 

dw = dz — adz, u = —4 r. (4.2) 

^ ^ l-de{z,z) ^ ' 

The Beltrami differential defines the new holomorphic structure - the deformed antiholo- 
morphic operator annihilating dw, while the antiholomorphic structure is kept unchanged 

= 9 + [id, du, = d. 

In addition, assume that /i vanishes in the marked points 

/i(z,z)U„ = 0. (4.3) 

In our construction we consider small deformations of the basic complex structure 
{z,z). It allows to replace ( [4.2|) by 

H = de{z,z). (4.4) 



Nevertheless, in some cases we will use the exact representation ( [4.2|) as well. 
4.2 Fibers 7^'. 

Let £^ be a principle stable G bundle over a Riemann curve „ . Assume that G is a 
complex simple Lie group . The phase space H' is constructed by the following data: 

i) the affine space {A\ of Lie(G)-valued connection on £. 
It has the following component description: 

a) connection {^4} , corresponding to the dw = dz component of A; 

b) The dual to the previous space the space {A} of dw components of connection A. A 
can have simple poles in the marked points. Moreover, assume that A/j, is a well defined 
function ; 

ii) cotangent bundles T*Ga = {{Pa,ga), Pa e Lie*(Ga), ga e Ga}, (a = 1, . . . , n) in the 
points (xi, . . . , Xn)- 

There is the canonical symplectic form on TZ' 



n 

j <5A,5A> +2'KiY,5 <pa,ga^5ga>, (4.5) 



where < , > denotes the Killing form on Lie(G). 

According with Sect. 2 we can consider the elements ga as trivializations of fibers in 
the marked points and Pa as the residues of holomorphic connections in the same points 

ga e Isom{Va, V), {V ~ Lie(G)), pa = res\oc,A. 
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Thus, the symplectic form ( [4 .51 ) is the generahzation of ( f^.lD on the singular curves. 
4.3. Extended phase space V . 

According to the general prescription the bundle V over M' with TZ' as the fibers plays 
role of the extended phase space. Consider the degenerate form on V' 



UJ = OJq 



- - / < 5 A, A > 511. (4.6) 



Thus, we deal with the infinite set of Hamiltonians < A, A > {z,z), parametrized by 
points of Sg „ and corresponding set of times fi{z,z). We apply the formalism presented 
in the previous section to these systems. 

4.4. Equations of motion 



They take the form (see (|373| ) , (|4^ )) 



— {z,z)=0, K—{z,z)=A{z,z), 

OjJL OjjL 

|5^ = 0. ^ = 0, (4.7) 

0/i on 



We can introduce the modified connection 



A' = A- -fiA. {U 

K 



In its terms ([4.6|) take the canonical form: 

» n 

uj= / <5A,5A' > + Y,5 <pa. g~^5ga > . (4.9) 



a=l 

and the equations of motion ( |4.7|) become trivial 

^^(.,.)=0, ^(.,.-) = 0. (4.10) 

In other words, as we discussed in Sect. 2, if we forget about the complex structures of 
curves the bundle become trivial P' ~ 7^' x AT'. 

4.5. Symmetries 

The form u ( |4.6|) (or ( [4.9|) ) is invariant with respect to the action of the group Qq of 
diffeomorphisms of „, which are trivial in vicinities Ua of marked points: 

go = {z^ N{z, z), z ^ N{z, z), ), N{z, z) = z + o{\z - xj), z e WJ. (4.11) 

In particular, its action on the Beltrami differentials takes the form of the Mobius trans- 
form 



/^-^^ (4-12) 



dz _j_ dz 
dN f^dN 
dz I dz 
ON r'dN 
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Another infinite gauge symmetry of tlie form ( [4 .61 ) (or ( |4.9| )) is the group 

^i = {/(^,^)gC°°(S„G)}, 
that acts on the dynamical fields as 

A + Kd ^ f{A + Kd)f-\ A + d + fid ^ f{A + d + i2d)f-\ 

{A' + d^f{A' + d)r'), (4.13) 

Pa ^ faPafa^, 9a ^ 9afa^, {fa = hm f{z, z)) , ^ fi. 

In other words, the gauge action of Qi does not touch the base Af' and transforms only 
the fibers 71'. The whole gauge group is the semidirect product 

Gi0Go- (4.14) 

There is an additional finite-dimensional symmetry group Q2 which commutes with 
( [4.14| ). It acts only on the singular curves in the fibers at the marked points. It is remnant 



of Qi on the desingular curves (see Sect. 2). 

Pa^Pa, Qa^haQa, ha E Ga, {a = 1, . . . , u) . (4.15) 

This action commutes with ([4.14|) . 

4.6. Symplectic reduction with respect to Qi 

The infinitesimal action of Gi ( [4.13| ) generates the vector field ei G A^^\T,g n,^ie(G)) 



Let 



je^A = -Kdti + [ei, A], 

je^Pa = [ei{Za,Za),Pa], 
je^ga = -gaei{z,Za), 
Jei/^ = 0. 

F^^^, = dA- kOA' + [A, A] = {d + dij)A - kOA + [A. A]. 



Note that the action of the operator {d + dn) on the connection A is correctly defined. 
The corresponding moment map 

Ji-.V^ Lie*(6?i). 

takes the form 

n 

Ji = -FaA^, z) + 2m 5\xl)pa. 

a=l 
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Let J7i = 0. In other words, the moment constraints equation takes the form 



F^^^,{z,z) = 2mJ2S'{x'i)pa. (4.16) 

a=l 

It means that we deal with the flat connection everywhere on Eg^„ except the marked 
points. The holonomies of {A, A) around the marked points are conjugated to exp 2iTipa. 
Let (L, L) be the gauge transformed connections 

A = fLf-' + f{d + fid)f-\ (4.17) 

A = fLf-' + Kfdf-\ (4.18) 

Then ([4.16|) takes the form 

n 

{d + dfx)L - KdL + [L, L] = 2m ^ 5^(a;°)p„. (4.19) 

a=l 

The gauge fixing allows to choose A in a such way that dL = 0. In fact, the anti- 
holomorphity of f~^{d + nd)f + f~^Af amounts to the classical equations of motion for 
the Wess-Zumino-Witten functional Swzwif, A) in the external field A, which does has 
extremal points. Then instead of ( [4.19| ) we have 

{B + dfi)L + [L, L] = J2 S\xa)pa. (4.20) 

a=l 

We can rewrite it as 

n 

d^L+[L,L] = 2mJ2S\xa)pa. (4.21) 

a=l 

Anyway, by choosing L we fix somehow the gauge in generic case. The last form of the 



moment constraint ( |4.21| ) coincides with that for the Hitchin systems ^ , which allows to 
apply the known solutions. 

4.7. Symplectic reduction with respect to Q2 

The gauge transforms ha G Ga in the points Xa acts on T*Ga- In the case of punctures 
it allows to fix pa on some coadjoint orbit pa = gaP^a^da^ and obtain the symplectic 
quotient Oa = T*Ga/ /Ga- In fact, the moment corresponding to this action is 

fJ-a = 9aPa9a^ £ Lie*(G). 

Let fia = Ja be some fixed point in Lie*(G'). Then the gauge fixing allows to choose ga up 
to the stabilizator of Ja and 

Pa = 9a^Ja9a G Oa- 

Thus, in ( [4. 19] ) or ( [4.20[ ) pa are elements of Oa- The symplectic form on Oa keeps the 
same form as on T*Ga- 

Consider the double points case (the nodal singularities). Let xi and X2 be preimages 
of the nodal point y under the normalization map. Then the symplectic form on the 
normalization of the singular curve 

uj = 6 < Pi,9i^Sgi > +5 < p2, 92^^92 > 



16 



generates the moment 
Put yU = 0. Then 



= 9iPi9i ^ + 92P292 ^■ 



Pi = -92P292 ^ 92 = 9i ^92- 



The pair {p2,92) is an arbitrary element of T*G. Therefore, in this case 
T*Gi®T*G2//G = T*G. 

In what follows we will concentrate on the case of curves with punctures (without 
double points). Let Ig^n be the equivalence classes of the connections {A, A) with respect 
to the gauge action ( [4.17|) , ([4.18|) - the moduli space of stable flat G bundles over Sg,^. It 
is a smooth finite dimensional space. If we fix the conjugacy classes of holonomies (L, L) 
around marked points „ becomes a symplectic manifold. It is extended here by the 
coadjoint orbits Oa in the marked points Xa, {a = 1, . . . ,n) in the consistent way (see 
( [4.19| )). Fixing the gauge we come to the symplectic quotient 

n 

^ C Ig,n X n C'a, 
a=l 

-1/ 



n = n'/fiGi © 02) = JT'(o)/(6?i © 02). 



It has dimension 



dim(ELia//G), 9 = 0, 

dim(7^) = <( 2rankG' + dim(ELia//i^), 9 = 1 
{2g - 2) dim G + dim(ELi a) ^7 > 2, 

where H is the Cartan subgroup H G G and Oa/ /G and Oa//H are the symplectic 
quotients of the symplectic spaces Oa under the actions of the automorphisms of the 
bundles in the rational and the elliptic cases correspondingly. The connections (L, L) in 
addition to p = (pi, . . . ,p„) depend on a finite even number of free parameters 2r 

(v, u), V = (t;i, . . . , W^), U = (Ui, ...,Ur). 

(7 = 0, 

rankG, 9 = ^1 

((7-l)dimG, g>2. 



7^ is a symplectic manifold with the non degenerate symplectic form which is the reduction 

of 

» n 
a=l 

On this stage we come to the bundle V" with the finite-dimensional fibers TZ over the 
infinite-dimensional base M' with the symplectic form 



» n 

/ <5L,5L> +2TnY.5 <pa,ga5ga' >= (4.22) 



IxJ 



I <5L,5L>+2m^5<pa,9a^5ga>-K,[ <L,6L>6fi. (4.23) 
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4.8. Factorization with respect to the diffeomorphisms Qo 

We can utilize invariance of u with respect to Qo and reduce A/"' to the finite-dimensional 
space Af, which is isomorphic to the moduli space Mg,n- Let cq be a vector field gen- 
erated by the diffeomorphisms ( ^.11| ). Consider the action of the Lie derivative £^0 
A=iA,A) 



C,,A = dj,,A + j,,dA = dAUeoA) + j,„F4, (F4 = dA+-[A,A]). 



Thus for the flat connections F4 = the action of diffeomorphisms Qq on the connection 
fields is generated by the gauge transforms jeo-^ ^ Gi- But we already have performed the 
symplectic reduction with respect to Qi. Therefore, j^,, belongs to the kernel of u (|4.23|) 
and we can push it down on the factor space V'/Qq. Since Qq acts only on M', it can 
be done by flxing the dependence of fj, on the coordinates in the Teichmiiller space Tg n. 
According to ( [4 .41) represent ^ as 

I 

^J' = ^^J's, fis = ts^J^^s, I = dimTg^n, ^J'°s = dris. (4.24) 

The Beltrami differential ( [4.24|) deflnes the tangent vector 

t = {ti, . . . ,ti), 

to the Teichmiiller space 7^^„ at the flxed point of Tg n- 

We specify the dependence of /i on the positions of the marked points in the following 
way. Let D Ua be two vicinities of the marked point Xa such that fl = for a 7^ b. 
Let Xa{z, z) be a smooth function 

Introduce times related to the positions of the marked points ta = Xa — x^^. Then 

fla = dna{z,z), na{z,z) = {1 + Ca{z - Xa))Xa{z,z) (4.26) 

In other words na{z, z) deflnes a local vector fleld deforming the complex coordinates only 

W = Z- taUaiz, z). 

The action of Qq on the phase space V" reduces the inflnite-dimensional component 
M' to Tg^n- After the reduction we come to the bundle with base 7^^„. Substituting 

5/x = ^/x°5t„ (9, = 9,J. (4.27) 

s=l 

in (|4.23|) we obtain 



1 ' 



t^o(v, u, p, t) ^ 5i^s(v, u, p, t)5ts, (4.28) 
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where ujq is defined by ( [4.22D , and Hg are the Hamihonians 

Hs= I <L,L> dsii (4.29) 
Jt, 

In fact, we still have a remnant discrete symmetry, since uj is invariant under the 
mapping class group TioiGo). Eventually, we come to the moduli space M.g,n = T^g,nl'^o{Go)- 

Summarizing, we have defined the extended phase space (the bundle V) as the result 
of the symplectic reduction with respect to the Qi®Q2 action and subsequent factorization 
under the Qq action. We can write symbolically 

v = {v"//gi®g2)/Go. 

It is endowed with the symplectic form ( [4.28| ). 

4.9. The hierarchies of the isomonodromic deformations (HID) 

The equations of motion can be extracted from the symplectic form ([4.2^), as it 



was described in Section 3 (see ( p.3|) ). They will be referred as the hierarchies of the 
isomonodromic deformations (HID). This notion will be justify later. In terms of the 
local coordinates (13.3) takes the form 



K^s^r = {H„\}^^, ndsU = {Hs,u}^^, ndsP = {Hs,p}ujo (4.30) 

The Poisson bracket {•, -joj,, is the inverse tensor to uo. We also has the Whitham hier- 
archy ( |3.4| ) accompanying ( |4.3(]| ). It follows from ( |4.22| ) that the Hamiltonians Hs ( |4.29| ) 
commute. 

{Hr, Hs}^g = 0. 

Therefore, 

dsHr - drHs = 0, (4.31) 

and there exists the one form on Aig^n defining the tau function of the hierarchy of 
isomonodromic deformations 

6\ogr = --J2HJts. (4.32) 

The following three statements are valid for the hierarchy of isomonodromic deforma- 
tions p:30| ): 



Proposition 4.1 The flatness condition ( \4-l!^ and HID (W^^W o,re equivalent to the 
consistent system of linear equations 

{Kd + L)^ = 0, (4.33) 

{Kds + Ms)^ = 0, {s = l,...J = dimMg,n) (4.34) 

{B + fid + L)^ = 0, (4.35) 
where Mg is a solution to the linear equation 

dioMs - [M„ L] = KdsL - (4.36) 
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Proposition 4.2 The linear conditions l \4-!i^4) pi"ovide the isomonodromic deformations 



of the linear system ^.33[ ), (2111/' with respect to change the "times" on M.g,n- 



Therefore, HID ( [4.30|) are the monodromy preserving conditions for the linear system 



The presence of derivative with respect to the spectral parameter w E Sg .„ in the 
hnear equation ( [4.33|) is a distinguish feature of the isomonodromy preserving equations. 
It plagues the application of the inverse scattering method to this type of systems. The 
later means the solving the Riemann-Hilbert problem which amounts the reconstructing 
the pair (L, L) from the monodromy data. In [35| this technique was applied to calculate 
asymptotics of solutions for bundles over rational curves. Nevertheless, in general case we 
have in some sense the explicit form of solutions: 

Proposition 4.3 (The projection method.) The solution of the Cauchy problem of 
for the initial data v", u°, p*^ at the time t = t^ is defined in terms of the elements 
L^,L^ as the gauge transform 



Lit) = f-\L\fi{t) - /i(t°)) + L')f + f'\d + /i(t)5)/, 

L{t) = f~\d + L')f, P(t) = r^(p°)/, 
where f = f{z,z) is a smooth G -valued functions on T^g.^ fixing the gauge. 



(4.37) 
(4.38) 



It means that solutions of HID are gauge transformations of free motion in the upstairs 
system. Equations ( 4.37|) , ( [4.38 ) look like the dressing transform of the free motion. To 
find solutions one should know the gauge transform / from the upstairs system to a 
fixed gauge. For example, in case of genus zero we consider the holomorphic solutions \E'. 
Therefore, / should kill the L operator in ( [4.371 ) (see below). 

Proois. 

To prove first statement represent A as ( 4.18 ). The first equation in ( [4.7| ) dsA = means 
that 

KdsL - ndMs + [M„ L] = 0. (4.39) 

Then ( [4.39| ) is the consistency condition for the linear system ( [1.33[) , ( [1.34[) . On this 
stage Ms is defined as Mg = nf'^dsf. To find the linear equation ( [4.36[) , defining Mg, it 
is necessary to substitute the gauge transformed form of A ( 4.17 ) together with ( [4.1^ ) in 
the equation of motion in the form ndgA = Adg^i- This equality is the same as 

K{^s^l^ + dsL) - d^Ms + [Ms, L] = dsfi{Kd + l). 

The later equation coincides with ( [4.36| ). Eventually, 



the compatibility ([CTI ), ( ^ ) 
the compatibility (1031) , ( F35| ) 
the compatibility ( |4.34[ ), 



the Lax equations ([4.39[ 



the flatness ([4T9[) 



This concludes the first statement. 
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To prove the second statement note that ( [4.3^j| ),( ^755D are equivalent to 

(Kd + A)"^^ = 0, 

{d + i2d + A)^f = 0, (^/ = /-i^). 

Due to the equations of motion {dgA = 0, dgA = ^Ads^) the monodromies of this system 
are independent on moduh, i.e. dg"^-^ = 0. The monodromies of the reduced system 
( [4.33| ),( ^l35D are conjugate to the monodromies of the later one. Thus, we come to the 



second statement. 

To derive the expressions for dynamical variables in the projection method we lift the 
initial data L°, L^,Pi, ■ ■ ■ ,Pn from the reduced phase space 71 in the point t° to 71' by the 
trivial gauge transform. Due to the equations of motion ( [4.7|) , the evolution in 7Z is trivial 



A{t) = L\ A = KL\^i{t) - ^i{t^)) + (L°) 

and can be push back on the reduced phase space 7Z by the gauge transform to the fixed 
gauge. This procedure is reflected in projection method formulae ([4.37|),(^4.38|). 



5 Remarks about the Hitchin systems and the KZB 
equations 

5.1. Scaling limit. 

Consider our system in the limit k — > 0. We will prove that in this limit we come to 
the Hitchin systems, which are living on the cotangent bundles to the moduli space of 
holomorphic G-bundles over The value k = is called critical and looks singular 

(see (^^), (|4.28| )). To get around we rescale the times 

t = T + Kt^, (5.1) 

where t''^ are the fast (Hitchin) times and T are the slow times. Therefore, 

s 

After this rescaling the forms ( [4.6D , (|4.28|) become regular. The rescaling procedure means 
that we blow up a vicinity of the fixed point (O) in Mg^n and the whole dynamic of the 
Hitchin systems is developed in this vicinity 0. For example, we have instead ( [4.5| ) and 

(ED 



UJ 



I < 5 A, 5 A > +2m 5 < g;^6ga > - E / <5A,A> dsfi^t"". (5.2) 

•^^ a=l s "'S 



If K = the connection A behaves as the one-form A G f2^^''^^(Sg^„, Lie{G)) (see ( [4.13| )). It 
is the so called the Higges field in terms of 0. An important point is that the Hamiltonians 
now become the times independent. The form ( |5.2| ) is the starting point in the derivation 

^We are grateful to A.Losev for elucidating this point. 
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of the Hitchin systems via the symplectic reduction |^. Essentially, it is the same 
procedure as described above. Namely, we obtain the same moment constraint ( [4 .201) and 



the same gauge fixing ( [4.17|) . But now we are sitting in a fixed point /i(t°) of the moduli 
space A^g,n and don't need the factorization under the action of the diffeomorphisms and, 
thereby, do not worry about the modular properties of solutions (L, L) of the moment 
constraint ( [4. 201 ). This only difference between the solutions {L,L) and the quadratic 



Hamiltonians Hg in the Hitchin systems and in HID. The symplectic reduction allows to 
identify the phase space TZ with the cotangent bundle to the moduli of holomorphic stable 
bundles. 

Propositions 4.1 and 4.3 are valid for the Hitchin systems in a slightly modified form. 

Proposition 5.1 There exists the consistent system of linear equations 

(A + L)^ = 0, AeC (5.3) 

{ds + M,)^ = 0, 5, = (s = 1, . . . , / = dim Mg,n) (5.4) 

((9 + L)^ = 0, 8 = d,, (5.5) 
where Mg is a solution to the linear equation 

BM, - [Ms, L] = dsL - L/. (5.6) 

Here we have 

The compatibility (|5.3|) , (|5.4| ) < — > the Lax equation Ot^L = [L,Ms], 
the compatibility ( |5.3| ), ( |5.5| ) < — > the moment equation ( [4.20| ), 



the compatibility (|]D, (U) < — > (|5]|). 



To derive these equations from the general case ( [4. 331 ), ( f4.34| ), and ( [4.35| ) we use the 
WKB approximation 

\1/ = <l>exp — , 

K 

where $ is a group valued function and 

S = 5(t, w) = 5°(T, w^) + kS\ = z-Y, TsUsiz, z), 

s 

S' = Y.tf{-^-ns{zrz)d)S'. (5.7) 

s C>T, 



In the first order in k the equations (^:^ , (|0^) , (|0^) just gives (^,(0),(^ if 



95" = A, (5.^ 
(9 8 



dwo ' dt^ 
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Therefore, as soon as iS° satisfies these equation, \1/ provides solutions to the hnearization 
of Hitchin systems. The detaile analyses of the perturbation in the rational case was 



undertaken in ^6 



Equation ( [5.3| ) allows to introduce the fixed spectral curve 

C: det(A-/rf+L) = 0, CeP(r*S©l) 

where A is a coordinate in the cotangent space. The Hitchin phase space TZ has the 
"spectral" description. It is the bundle 



c 



over the moduli A^c" of spectral curves with abelian varieties as generic fibers. The map 
TT acts from the pair [A, A) ~ (L, L), (L now is the Higgs field L G fi(i'°)(Sg,„, Lie*(G'))) 
to the set of coefficients of the characteristic polynomial det(A ■ Id + L). The one-form 
6 = Xdw^ being integrated over corresponding cycles in C gives rise to the action variables. 
The angle variables can be also extracted from the spectral curve. All together defines 
the symplectic structure on the Hitchin phase space in the spectral picture. This original 
Hitchin construction is working in the singular well 0. The symplectic structure 

of this type connected with hyperelliptic curves was introduced originally in the soliton 
theory by Novikov and Veselov 1^^. In terms of 4d gauge theories 6 is the Seiberg-Witten 



differential. It follows from ( |5.8| ) that 6 = dS^ and along with ( p.7| ) it defines the first 
order approximation to the solutions of the linear form (|4. 33|) , (|4 . 34| ) , ( ^735]) of HID. 



It is possible to define the dynamic of the spectral curve beyond the crtical level as it 
was done for the Painleve equations in It follows from the Lax representation (|4.39|) 



that 

ds det(A ■Id + L)= tr(9^M,(A ■ Id + L)-\ 

Note, that it defines the motion of C only within the subset M.g^n C M.c- 

When L and thereby M can be find explicitly the simplified form of ([4.20|) allows 
to apply the inverse scattering method to find solutions of the Hitchin hierarchy as it 
was done for SL(A^, C) holomorphic bundles over Si i ||3^, corresponding to the eUiptic 
Calogero system. We present the alternative way to describe the solutions: 

Proposition 5.2 (The projection method.) 

Lit) = f-'L'f, p^it) = r\pl)f 



The degenerate version of these expressions was known for a long time |40 



5.2. About KZB 

The Hitchin systems are the classical limit of the KZB equations on the critical level 
11^, 10|. The later has the form of the Schrodinger equations, which is the result of 



geometric quantization of the moduli of flat G bundles [0, The conformal blocks of 
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the WZW theory on with vertex operators in marked points are the ground state 
wave functions 

HsF = 0, (s = l,...,/). 

The classical limit means that one replaces operators on their symbols and generators of 
finite-dimensional representations in the vertex operators by the corresponding elements 
of coadjoint orbits. 

Generically, for the quantum level k"^"""* 7^ the KZB equations can be written in the 



form of the nonstationar Schrodinger equations |T^, |T0 



0. 



To pass to the classical limit in this equation we replace the conformal block by its 
quasiclassical expression 

F = exp — , 

where S is the classical action {S = logr ( [4.32D ) and renormalize 

l^quant 



h 

The classical limit ^ — > 0, k''"'*"* leads the Hamilton- Jacobi equations for S, which 
are equivalent to the HID ( ^4.30| ). 

Summarizing, we arrange these quantum and classical systems in the commutative 
diagram. The vertical arrows denote to the classical limit and mean the simultaneous 
rescaling of the quantum level, while the limit k''"""* — ^ (k — >■ 0) on the horizontal 
arrows includes also the rescaling the moduli of complex structures. The examples in the 
bottom of the diagram will be considered in next sections. 



KZB eqs., {K,Mg^n,G) 
(a = 1, . . ., dim Mg^n) 



^quant 



KZB eqs. on the critical level, 
(a = 1, . . . ,dimA^g_„) 



Hierarchies of Isomonodromic 
deformations on Aig^n 



Schlesinger eqs. 
Painleve type eqs. 
Elliptic Schlesinger eqs. 



EXAMPLES 



Hitchin systems 



Classical Gaudin eqs. 

Calogero eqs. 
Elliptic Gaudin eqs. 



6 Genus zero - Schlesinger equations 

6.1. Derivation of equations 
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Consider CP^ with n punctures (xi, . . . ,x„|xa 7^ Xb)- The Beltrami differential fi is 



related only to the positions of marked points. Then from ( [4.26|) 

6fia = d{l + Cad{z - Xa)Xa{z, z))6ta, {Sta = SXa). 



(6.1) 



On CP^ the gauge transform ( ^.17| ) allows to choose A to be identically zero. After 
the gauge fixing 

A = fd^f-\ (6.2) 
A = fLf-' + Kfd^f-\ 
the moment equation takes the form 

n 

du,L = 2m^5'^{Xa)Pa- 



a=l 



It allows to find L 



Then we have from (|0| 



Pa 



\W-Xa 



(6.3) 



-5 I < L,L > 6a 
2 Jcpi 



< Pa,Pb > 



icpi (W - Xb) (W - Xa) 



J2 d{l + Cc{w - Xc))XcSXc = ^(5iJa,l + 5Hafi)5Xa, 



where 



and 



Ha,, = E 



< Pa,Pb > 



(6.4) 



H2,a = Ca< Pa, Pa > ■ 



Hi a are precisely the Schlesinger's Hamiltonians. On the symplectic quotient uo (|4.23|) 
takes the form 

n n 

UJ = "^6 < Pa, Qa^^Qa > '^{6Hb,l + 5Hbfi)5Xb. 

a=l ^ 6=1 

Note, that we still have a gauge freedom with respect to the coordinate independent 
G action, since this action does not change our gauge fixing (|6.2|) . The corresponding 
moment constraint means that the sum of residues of L vanishes: 



0=1 



(6.5) 



While H2,a are the Casimirs and lead to the trivial equations, the equation of motion 
for Hi a are the Schlesinger equations 

r. [Pa,Pb] , / , X 

i^ObPa = — - — , [a t b), 

Xa Xb 
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I^OaPa = - 2^ 



As by product we obtain by this procedure the corresponding hnear problem ([4.33|) , (|4.34|) 
and (|4^ ) with L = Q, Z = and 



W - Xa 



as a solution to ([4.36|) . The tau-function ( [4. 321) for the Schlesinger equations has the form 



(51ogr = - ^ < > '51og(xc - Xb). 

ci^b 

6.2. Solutions via the projection method 

We will find the dressing transform defining the evolution on the coadjoint orbits 

p,(t) = ri(2;,^,tK/(2;,^,t). (6.6) 

Recall that the times in the Schlesinger equations are ta = Xa — x^ and assume that t° = 0. 
We have from ( |6.3| ) 

n 

L(to = 0) = L° = E 



a=l ^ 

It follows from ( |6.2| ) that L = 0. Then the projection method ( [4.37D gives in this case 

r\z, z, t)LV(t)/(;2, z, t) + f-\z, z, t){d + fi{t)d)f{z, z, t) = 0. 

In other words, the gauge transform f{z,z,t) defining the evolution of solutions (the 
dressing transform) can be found from the equation 

[d + Yl tkdxkiz, z){d + L^)]f{z, z, t) = 0. (6.7) 

k 

We seek for smooth solutions to this equation assuming that the times tg = Xg — are 
small. To this end consider the perturbative series 

f{z,z,t) = id + '^tkak + '^tjtkajk + Utjtkaijk + ■ ■ . , (6.8) 

k j'^k i^j^k 

where ak,ajk,aijk, ■ ■ ■ are smooth maps CP^ Universal enveloping algebra(G'). Then 
( |6.7D leads to the system 



1) dak = -L^dxk 

2) dajk = -{d + LP){dxjak + dxka-j) 
S)daijk = -{d + L^)dx[jCijk] 
4)..., 

where [...] means the symmetrization. All equations have the same structure - their solu- 
tions depends only on the previous step. Since (0, l)-forms on CP^ are exact the equations 
can be integrated and solutions are found step by step. 
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In the first order one has 

Pa 



ak{z,z) = -J2 ^Xkiz,z), {xt{z,z) = Xk{z,z) - Xk{Xa,Xa)). (6.9) 



„ . z — X 

a=l 



Note that ak{z,z) is a nonsingular function on CP^ due to the definition of Xki^,^)- 
Consider now the second order approximation. 

da,,{z, z) = Y: d{-^)dxiA^, z)xt]{z, z) + Y: -^BxiA^, z)dxt]{z, z) + 

a=l ^ X^ ^^-^ Z 

n n _ 

+ E-^E-^dx,i^rz)Aiz,z). 



a=l ^ -^a h=l ^ 



The result of integration is 

. -^ o/M^llM^ x']{z,z)xt{z,z) 
ajkiz,z) = LPa( \_^o (^_^o)2 )+ (6-10) 

n ( 0\2 n 

E r^r^2xA^^ ^)xt{z, z)-Y. -^ix-iz, z)xt{z, z) - x';{xi)xi{xi)). 

a=l ^a) a^b ^ -^b 

Here ipj{z, z) is defined as the resuh of integration 

dr,{zrz) = 2xjiz,z)x';iz,z), {i;]{xa) = 0). (6.11) 

It provides the absence of poles in z = Xa in the first term in the right hand side of ( |6.10|) . 
Due to the subtraction of Xj{Xb)Xk{^b) the last term we kill the pole of ajk{z,z) in 
z = Xb- Thus, in second order approximation we obtain the regular solution as well. It is 



defined almost explicitly up to the integration (6.11 



Therefore, we have defined the dressing transformation (|6.8|) , (|6.9| ),( |6.10|) up to the 
third order of the initial data ( [6.61 ). The calculations of the higher order corrections 
reproduce the same procedure as on the second order and we can repeat them step by 
step. 



7 Genus one- elliptic Schlesinger, Painleve VI... 



The genus one case is still feasible to write down the explicit formulae for the Hamiltonians 

and the equations of motion. 

7.1. Deformations of elliptic curves 

In addition to the moduli coming from the positions of the marked points there is 
the elliptic module r, Imr > of curves As in (|4.24]) , (|4.26| ) we take the Beltrami 

differential in the form 



a=l 

where na{z, z) is the same as in ( [4.26| ) and 

n 

{Z- z){l~Y.Xa{z,z)). 
a=l 



(7.1) 
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We replace 

U 

tr ^ —, tr = T - To, p = Tq - Tq. 
P 

Here tq defines the reference complex structure on the curve 

Tq = {0 < X < 1, < y < 1, z = X + Toy, z = x + foy}. 
For small tr from ([7. 1|) 



fir = Prd{z - Z){1 -Y,X'aiz,z)), {^r = —) , (7.2) 

a=l T-TO 

or 

tr - _ 

fir = —d{z - Z){1 -Y.Xa{z,z)). (7.3) 
P a=l 

As we assumed from the very beginning, fir vanishes in the marked points. It describes not 
only a small vicinity of Tq in A^i i, but also the whole Teichmiiller space as well. In terms 
of fir the Teichmiiller space is the unity disk \fir\ < 1 (see ( |4.12| ) for the transformation 
law). In terms of r it is the upper half plane Imr > 0. On the other hand, the times 
t = {tr, ti = xi — x\, . . . ,tn = Xn — x°) define deformations 

^0 ("^l' • • • 5 -^n) Tr (3^1) • . . ) ^n) 

in a vicinity of Tq in M.i^i. Eventually, for small tr and ta we have from ( [7.3|) 



5fl = 5fird{z - Z){1-Y1 Xaiz, Z)) + ^ Stad{l + Ca{z - Xa))Xa{z, z)) (7.4) 



a=l a=l 



7.2. Flat bundles on a family of elliptic curves 

Note first, that A as a (9-connection determines a holomorphic G bundle S over T^. 
For stable bundles A can be gauge transformed by ( |4.17| ) to the Cartan {z, z)-independent 
form L 

A = fLr' + f{d+fid)f-\ 

L ^Ti — Cartan subalgebra of Lie(G'). 
Therefore, a stable bundle £ is decomposed into the direct sum of line bundles 

S = ®UiCk, r = rank(G). 

The set of gauge equivalent connections represented by {T} can be identified with the r 
power of the Jacobian of T^, factorized by the action of the Weyl group W of G. Put 

I = 27r^i^u, vleH, (1^ = ^^). (7.5) 

P P r -Tq 

It means that 

J^^ Ldw = u. (7.6) 



28 



Let L = d log 0. Then the integral 

rP 



/ Ldw = / log (f)dw. 
Jt3 JPn 



'T? JPo 

defines the Abel map P E ^ u. We will come again to this point later. 

The fiatness condition (the moment constraints ( [4.21| )) for the gauge transformed 
connections (L, L) takes the form 

n 

d^L + [L, L] = 27n ^ 5\xl )pa. (7.7) 

a=l 

Let R = {a} be the root system of of Lie(G) = Q and 

Q = 1-L (BaeR Ga 

be the root decomposition. Impose the vanishing of the residues in ( |7.7| ) 

n 

EPa|H = 0, (7.8) 

a=l 

where PaW is the Cartan component of Pa and we have identified Q with its dual space 
Q* . This condition is similar to ( |6.5|) and has the sense of the moment constraints for the 
remnant gauge action (see below 7.3). 

We will parametrized the set of solutions of ( [r.7|) by two elements v, u G 7i. Let Ei{w) 
be the Eisenstein function of module r ( [A. 2] ) and {pa)'H aiid {pa)a are the Cartan and the 
root component of Pa € Ca, 

Lemma 7.1 Solutions of the moment constraint equation ( \7.T{ ) have the form 

L = P + X, Pen, X=Y.X^. (7.9) 



V U 



P = 2m{-—^ -K- + J2{pa)HEiiw - xj), (7.10) 

^ pT P a=l 

n 

X^ = Y.Xl (7.11) 

a=l 

X^ = — exp 2Tn{ a{u)}(f){a{u),w - Xa), 

I- fir r - To 



where (j){u,w) is defined in ( \A.4) - 



Proof. Consider first the Cartan component. Since L Eli. 

n 

dij,P = 2Tli^5'^{Xa)Pa- 



a=l 
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From ( [A.28D ,( [A.30D we obtain ( |7.1(J| ). The special choice of the constant part of P will be 
explained later. Note, that v G ?i is a new parameter. 



For the root components ( [7.7|) takes the form 

.1 - /ir 



P 



a=l 



Comparing it with ( |A.31| ) and its solution ( |A.32| ) we come to ( [7.11[ ).D 
Therefore we have found the flat connections L = L(v, u), L = L[u). 

7.3. Symmetries. 

The remnant gauge transforms do not change the gauge fixing and thereby preserve 
the chosen Cartan subalgebra Ti G G. These transformations are generated by the Weyl 
subgroup W oi G and elements f{w,w) e Map(T^, Cartan(G')). Let 11 be the system of 
simple roots, = {a"^ = j^;^}, is the dual root system, and m = Z^Qgn^^aCt"^ be the 
element from the dual root lattice ZK^ . Then the Cartan valued harmonics 



exp 27rz(m- 



w 



W TW 

— + n 

r 



r - To t-Tq 
generate the basis in the space of gauge transforms. They act as 

m — nr 



(7.12) 



L + 2m- 



T-To 

-m + nfo 



P ^ P + 2mK- 

r - To 



exp 



■[{ma + naTo){w - Xa) - (m^ + naT){w - Xa)]- 



T -To 

In terms of the new variables v and u they take especial simple form 



u 



u + m — nr. 



KIl. 



(7.13) 



(7.14) 



The whole discrete gauge symmetry is the semidirect product W of the Weyl group W 
and the lattice ZK^ © rZR^ . It is the Bernstein-Schvartsman complex crystallographic 
group The factor space Ti/W is the genuine space for the coordinates u, that we 
discussed above (see ( [7151 ) and ( [7.6| )). 

The transformations (|7.12|) according with (|4.13|) act also on pa ^ Oa- This action 
leads to the symplectic quotient Oa/ 1 H and generates the moment equation (|7^). 

The modular group PSL2(Z) s a subgroup of mapping class group for the Teichmiiller 
space Ti^n- We don't consider here the action of the permutation of the marked points 
on dynamical variables (v, u, p, r, x^). Due to ( [4.12| ) and ( |7.2| ) its action on r takes the 
standard form 

r^^r = ^^ 7ePSL2(Z). 

cr + d 

We summarize the action of the Bernstein-Schvartsmann group and the modular group 
on the dynamical variables: 
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W={s} 


Zi?^ © tZR^ 


PSL2(Z) 


V 


sv 


V + /tn 


v(cr + d) — Kcu 


u 


su 


u — m + nr 


u(cr + d)^^ 


{Pa)'H 


s{.Pa)H 


Pa 


Pa 


{Pa)a 


iPa) sa 


^irna,na){Pa)a 


{Pa)a 


T 


T 


T 


ar+b 

CT + d 








CT + d 



Here ip{ma,na) is defined by ( [7.13| ). 



7.4 Symplectic form. 

The set (v, u) G 7i, p = (pi, . . . ,p„) G ©"^i^a of dynamical variables along with the 
times t = {tT.,ti, . . . ,tn) describe local coordinates in the total space of the bundle V. 
According with the general prescription, we can define the hamiltonian system on this 
set. The main statement, formulated in terms of the theta-functions and the Eisenstein 
functions (see Appendix), takes the form 



Proposition 7.1 The symplectic form ui ( \i-2Sl ) on V 



IS 



1 " 1 " 1 

UO = (5V, H + E ^ < Pa. 9a^^9a > — (E ^^2,a + m,a)^t, " (7.15) 



47r2 

with the Hamiltonians 



a=l 



K 



a=l 



K 



H2,a = Ca< Pa, Pa >] 



Hi^a == '^{t—-^ l^-,Pa\H) + J2iP<^\'H,Pb\H)E: 

1 - P t^a 



l{Xa - Xb) + 



(7.16) 
(7.17) 



^Y^{Pa\a,Pb\~a)(t){a{vi),Xa - Xfo); 



V, V 



(7.18) 



E Y^i.PaU,Pa\~o)E2{a{u)) + ^ (p„ | Pfc | w) (^2 (a^a - X^) - E^Xa - Xfc))- 
^TI" a=l a a^b 

-^^J2iP<^\a,Pb\-a)4>{-<yiu),Xa - Xb){Ei{a{u)) - Ei{Xb - Xa + "(u))); 



(a(uj,x„ - Xb) 



d{a{M)+Xa-Xb)e'{id) 

9{a{u))9{xa - Xb) 
Proof. The form we have to calculate is ( ^4.23| ) 



UJ 



. _ n ^ 

/ < 5L, 5L > +2Txi S <pa, Qa^Sga > — / < L,SL > 5n 
7s ^1 K 7s 



with Sfi 



First, we gauge transform (L, L) by 

w — w 

8xp(2m 

a=l 



'(JJ — '^J 

f{w,w) = Y[exp{2TTi x'a{^,w)u) 



T-Tq 
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where we choose x'ai'^y'^) ^ such way that 

SUppXaiuJ,w) C SUppXa(w,w), (7.19) 

suppxK"^)""^) n svippdxaiuj , w) = 0, 

and Xaiwjw) related to the moduh curves (|7. 1|) . In fact, the first condition follows from 
the second. 

As we know the gauge transformations do not change uj, but do change relations 
between its summands. Instead of ( |7.5| ).( [rTU ) and ( |7.11| ) we obtain 

- 1 - u 

L = 2m ^ud^iw -w){l-Y^ x'ai^, w)), (7.20) 

P a=l 

V U " 

P = 2m{- - + Y.{pa)HEl{w - Xa))- 

1 - fir - ^ I 

-"^d^iw -w)^Xa^w,w), (7.21) 

P a=\ 
\ yj — yj " 

= ^exp{ — a(u)} V(pa)a0(tt(u),ti' - Xa), (7.22) 



Taking into account the explicit form of L ( |7.5| ) we obtain 

<SL,SL>=^^^^^ + SiST,St), (7.23) 
P 

where S{6t, 6t) is a sum of terms with a linear dependence on the "time" differentials. It 
is compensated by terms coming from 



47r2 

Let us calculate the Hamiltonians 



4- / < L,6L> 6fi. 



1 V u " 

< L^L >= {y—z + PalnEiiw - Xa)f- 

- _ ~ 'YiPa\a,Pb\-a)(p{-a{u),W - Xa)(p{a{u) , W - Xfe)- (7.24) 

pT a,b aeR 

^ ~ - + X!(Pa)w^l(w - Xa)), ^u)d^{w - w) X'ai'>^,'w). 

To take the integral over Tq we have to couple this two form with 6p (|7.4|) 

-7-1 / <L,L>6p = 
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-— / < L,L > 6fird{w - W){1 -Y,Xa{w,w)) + Y,^tadin{'^ + Ca{w - Xa))Xa{w,w)). 
^"^ a=l a=l 

Due to our choice of x'ai'^j'^) (|7-19|) , the last hue of (|7.24|) does not contribute in the 
integral. Therefore, we leave with the holomorphic double periodic part of two-form 
< L,L >. Using (|A. 10|) we rewrite ( [7.24D as 



47r2 



< L,L >-- 



U V u 

, r K— ) + 



I- fir p I- Hr P' 



V u 

- K-,Pa\H)El{w - Xa) + 



P 



'^{Pa\H,Pb\H){Ei{w - Xa) - Ei{w - Xfe))- 



. ~ J2 (PaU,Pf,|-a)0(-tt(u), w - Xa)0(a(u), w - 0;^)- 

^ Pr ajtb adR 



An 



-—- Y: E iE2{w - Xa) - E,{aiu)). 

^ pT a aGR 



Since L has only first order poles ( |7.21| ),( [7I2^ ), we expand it on the deformed torus 
according with (|A.34 ) 



An" 



< L,L>= (E H2,aE2{w - Xa) + Hi^aEi{w - Xa)) + h. 



(7.25) 



a=l 



Due to ([70^ ),( |O9| ) and (pOOl) 



-7-^ / < L,L > 6p = 

n . 

E( / {H2,aE2{w - Xa)du,{w - Xa)Xa{w,w) + Hi^aEl{w - Xa))du,Xa{w , w))Sta 
a=l T 

n 

+ hQ5prd{w - W)(l - E Xa{w, W)). 



a=l 



Taking into account and (|A.35| ),( p\.36| ) we find 

H2,a = res\^^ < L,L > {W - Xa) = Ca < Pa, Pa > +COnst., 

Hi,a = resU„ <L,L>= ( fm|) . 

The constant term in ( [r.25| ) ho is defined by ( |A.37|) . To find it we use (|A.10|) and (|A.38 
Then = h^drPT- After some algebra we obtain ( [7.18| ). □ 



7.5 Example 1. PVI^ _^ ^ 



4 ' 4 ' 4 '2 4 
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Consider SL(2, C) bundles over the family of Si^i. Then ([T.SD takes the form 



1 — T , 
L = 2ni diag(M, —u). 

P 



(7.26) 



In this case the position of the marked point is no long the module and we put xi = 0. 



We have from ( |7.2| ) 



W = Z [Z — z), w = z. 



du, = d + - — —d. 

r - To 

Since dim O = 2 the orbit degrees of freedom can be gauged away by the hamiltonian 
action of the diagonal group. We assume that 

p = z/[(l,lf 

Then we have from(|7:9|),([7lO|),(fm]) 



L 



x{u, w, w) 



27Ti( 





1-flT P' 



x{2u, w, w) 



2m{l-fi)^ 
The symplectic form ( |7.15| ) 

1 



X{-2U, W, W) 27Ti{-j^ + K^) 



exp 27Ti{{w — w)u —}(j){u, w) 



P 



„uj = (6v,6u) 6Hr6T, 

Stt^ k 



and 



Hr = — + U{u\t), U{u\t) = -{—yE2{2u\r) - — . 



2^' 



2m 



Then the equations of motion are 



dv 



du 

Z/2 d 



K 



Or {27riy du 



E2i2u\T). 



Due to ( |A. lOp we obtain 



z/2 d 



K 



dr"^ {2Txiy du 



p{2u\t) 



(7.27) 



(7.28) 
(7.29) 



which coincides with ( |1.12|) . This equation provides the isomonodromic deformation for 
the linear system ( [4.33D ,( |05| ) with L ( |7.27] ) and L ( [7. 261 ) with respect to change the 
module r. The Lax pair is given by L ( [7.27| ) and Mj. 



Mr 



y{2u,w,w) 
y{—2u,w,w) 
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where y{u^w^w) is defined by the equation (see ( [4. 361 ) 

{du, H —u)y{u, w, w) = -. —:^x{u, w, id). (7.30) 

r - To K[r - Toy 

Using the representation (|A.33|) for x{u,w,w) = g2{u,w,w) we find 

y{u, w, Hi) = . / — —dux{u, w, w). (7.31) 

The equivalence of the Lax equation 

drL - Kdy,M + [M, L] = 

to the equations of motion ( [7.28| ),( fr29| ) can be checked by substituting L ( |7.27| ) and M 
with y{u,w,wy (|7.31| ). Details of this procedure will be replenished for the SL(iV, C) 
bundles in the next example. 

The projection method determines solutions of (|7.28| ),( [7.29|) as a result of diagonal- 
ization of L ( |7.27| ) by the gauge transform on the deformed curve 

-diag(M(r), -M(r)) = 



r - To 



_ 27u(r-7o)/ v'-i^f x{2u',z,z)\ 1 _ 

+ f~\z,Z,T){d+^^d)f{z,Z,T). 

T-Tq 

On the critical level {k = 0) we come to the two-body elliptic Calogero system. 
7.6. Example 2. 

For flat G bundles over Si i we obtain PVI-type equations, related to arbitrary root 
systems. They are described by the system of second order differential equations for 
the u = {ui, . . . ,Ur), {r =rankG) variables. In addition, there are the orbit variables 
p G 0{G) satisfying the Euler top equations. Consider in detail the SL(A^, C) case with 
the most degenerate orbits O = T*CP^~^. They have dimension 2N — 2. The orbit 
variables can be gauged away by the diagonal gauge transforms and we are left with the 
coupling constant z/. We already have the L and L matrices (see ( |7.5| ) and ([7.9|)) 

L = P + X, 
P = 2m{ z K-), 

1- fir P 

_ W — W 

X = {xa} = (r — ro)^exp 27rz{ — «(u)}0(a(u), w). 

T-Tq 

Here a = ej — e^, a(u) = Uj — Uk, (j 7^ k), and 

u = diag(Mi, . . . , un), V = diag(t;i, ...,vn)- 
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The equations of motion take the form 



K, 



duj 



K 



dr 
1 



dr 

-- -a,f/(u|r) 



(7.32) 
(7.33) 



UHt) = J^^Y.P3,kPk,jE2iUj - Uk\T). 



On the critical level this Painleve type system degenerates into A^-body elliptic Calogero 
system. 

Let us check the consistency of the equations of motion with the Lax representation 



drL - Kd^M + [M, L] = 0. 
Put the M operator in the form 



(7.34) 



M 



-D + Y, 



where 



Y = {yiUj-Uk)} = {Va}, 

and D is a diagonal matrix. We have already found i/a ( |7.31| ) using ( (4.361 ). But the 
diagonal part D is not fixed by ( [4.36|) and should be found from the consistency of the 
Lax equation with the equations of motion. We take it in the same form as in the Calogero 
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D = diag((ii, ...^d^), dj = Y^ s{uj - Ut). 

First, we will prove few facts concerning the matrix elements of L and M. We will 
prove that they satisfy the following functional equation 



x(m, 2, z)y{v, z, z) — x{v, z, z)y{u, z, z) = (s(f ) — s{u))x{u + v, z, z). 
In particular, its solution s{u,w) is w independent 

siu) = —p(u) + const. 



(7.35) 



(7.36) 



The relation ( [7. 351 ) is the so called Calogero functional equation. Due to ( |7.31| ) we can 
put in it the derivatives of x instead of y. The exponential factor in the expression of 
X cancels and ( |7.35| ) take the form of the addition formula ( |A.26| ). Simultaneously, we 
obtain ( [7.36|) . We also need the following identity 



{dr - ^d^du H ^9„)x(m, z, z) = 0. 

2m T — To 



(7.37) 



It can be derived from the representation ( [A.33| ) for x{u, z, z) = g2{u, z, z). 

Consider now the Lax equation (|7.34 ). Separation the diagonal and nondiagonal terms 
leads to the system 

-^P + tid^D + Y^{yc,x.a-y~aXa) (7.38) 



36 



—Xa - t^duiVa - a{D)xa + [V', X]^ = 0. 
Due to (|7.36| ) and ( [7.32|) the first equation ( [7.38|) can be rewritten as 

2'Ki-^Vk = /^^) ^{x{Uj - Uk)x'{Uk - Uj) - x{Uk - Uj)x'{Uj - Uk)) 



(7.39) 



dT 



2mK 



But 



x{u)x{—u) 



Y^[x{u^ - Uk)x{Uk - Uj)]'. 



(m, w)(t>{—u, w) 



iE2{w)-E2{u)), 



(see ( |A.10|) ). Thereby we come to the i VM) . 

Now check the off-diagonal part ( |7.39|) . It is convenient to go back from 



w 



r - To _ _ _ 
z, w = z 



to {z, z) variables. Then ( [7.39|) takes the form 

-^x{Uj - Uk) - KdzViuj - Uk) - 2ni{- — 

UT i — fi-r 



Uj Uk , I N 

K )y[uj - Uk)- 



N 



N 



—xiu 



-Uk)[J2siuj 



Ui 



^s(Mfc - Ui)] + 



N 



+ yi'^i ~ Ui)x{ui - Uk) - y{ui - Uk)x{uj - Ui) = 0. 



i=l 



It follows from ([7:31|) , (^) , and (HTSTj) that 

Uk) - Kd^y{uj - Uk) - 2TTi{^ — ^ - K— — —)y{uj - Uk) = 0. 



-x^u, 



p 



I- fir 

On the other hand the equality 

x{uj - Uk){s{uj - Ui) - s{uk - Ui)) + y{uj - Ui)x{ui - Uk) - y{ui - Uk)x{uj - Ui) = 0. 

is the addition formula (|7.35|) . It concludes the proof of the equivalence of the Lax 
equation and the equations of motion. 

The new ingredients of this construction in compare with the original form are the 
dependence the L-matrix on the spectral parameter^ , the presence of derivative d^M in 
the Lax equation, and the replacement of the external time t on the modular parameter 
r. Nevertheless, the form of the Lax matrices is defined as for the Calogero system by 
the same functional equation ( [7.35| ). It turns out that its solutions satisfy additional 
differential equations ( [7.36| ), ( |7.37| ), which allows to apply them in the isomonodromic 

■^The dependence on the spectral parameter first introduced in |^ for tlie Calogero system in a slightly 
different form. 
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situation as well. Again, the solutions of the equations of motion ( |7.32| ), ( [7.33D are 
obtained by the diagonalization 

27rzl^diag(ui, ...,Un)= lirr\L{^,\ u", r°)/+ 
P 



Appendix A 



We summarize the main formulae for elliptic functions, borrowed mainly from |^ . We 
assume that q = exp27rir, and the curve is C/Z + rZ factor of C under the shifts 
generated by (1,t). 

The basic element is the theta function: 

e{z\T) = q-^Yl (-l)"e"*W'^+i)-+2„.) _ 

neZ 



-2i'KZ\ 

L - q )yL - q o )yL - q o J 
n=l 

The Eisenstein functions. 



Ei{z\T) = dAoge{z\T), ^i(^|r) -- + ..., (A.2) 

z 

E2{z\r) = -d,E^{z\T) = dl loge(z|r), E2{z\r) ~ ^ + . . . . (A.3) 



z^ 



The next important function is 



(l>{u,z)= \ ' ; ' . A.4 
6[u)6[z) 

It has a pole at 2; = and 

res\^=Q(t){u, z) = I. (A. 5) 

Relations to the Weierstrass functions. 

C{z\T)=E^{z\T)+2ri,iT)z, (A.6) 

piz\T)=E2{z\T)-2r]^{r), (A.7) 

where 

Viir) = al)= (A.8) 
3_ ^ ^ 1 2Ar]'{T) 



where 

r^(r) = g2^[](l-g«). 



n>0 
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the Dedekind function. 

a(u + z) 
= exp{-2r]iuz)———. 

a[u)a[z) 

(f){u, z)(f){-u, z) = p{z) - p{u) = E2{z) - E2{u). 
Series representations 



Mz\r)--2m{- + ^^^) 



2 ^n^2mz 

n<o y ^ n>o y ^ ^ 

£;2(^|r) = -47r2 ^ 



y ^ ) 



g— 27rmz 



neZ 



Parity. 



e{-z) = -e(^) 

E^{-z) ^ -E^{z) 

E^i-z) = £;2(^) 

z) — (l){z, u) — —(f){—u, —z) 

Behaviour on the lattice 

e{z + 1) = -e{z), e{z + r) = -q-^-^^^'^eiz), 

Ei{z + 1) = Ei{z), Ei{z + r) = Ei{z) - 2ni, 
E2iz + l)^E2iz), E2iz + T)^E2iz), 

(j){u + 1, = z), (j){u + T,z) ^ e~'^'"'-''(j){u, z). 
Modular properties 

"(^TTrf l^7Td> = ''"('^ + axp(— )«(.-|r), (e» 

„ , z ,aT + b. , ,v ^ / , X 

£^1 ^ ^ = [cT^d)E^{z\T) + 2mz. 

£^2(^-^1^^^) = {CT + dfE2{z\T) + 2ni{cT + d) 
Addition formula 

(t){u, z)dy(f){v, z) - (f){v, z)du(f){u, z) = {E2{v) - E2{u))(l){u + v, 
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or 

z)d^(f){v, z) - (f){v, z)du4>{u, z) = {p{v) - p{u))(f){u + v,z). (A.26) 

The proof of (HKM ) is based on ([AlSD , (|AT7|) , and (^^T]). 

In fact, (f){u, z) satisfies more general relation which follows from the Fay three-section 
formula 

Zi)(t>{u2, Z2) - (f){Ui + U2, Zi)(f){u2, Z2 - Zi) - + U2, Z2)(f){ui, Zi - Z2) = (A.27) 

Green functions 

The Green functions are (1,0) forms on T^. The first gi{z) is defined by the equation 

dgi{z) = 2mJ2PaS\xa), (Ep- = 0), (A.28) 

a a 

where 

S^{Xa) = E fm,n{^ - - Xa), (A.29) 
m,n£Z 

27rz 

fmA^, z) = exp {m{z -z) + ni^rz - fz)}, {p = T-f), 

P 

9i{z) = "^PaEiiz - Xa) + const. (A.30) 

a 

For the equation 

(d+ — u)g2{u,z) = 2mS\0), (A.31) 
P 

the Green function is 

g2iu,z) = -c'''''^^<l>iu,z). (A.32) 
P 

g2{u,z)=p E -^^^^^^ (A.33) 

m,nGZ 

Expansion of elliptic functions. 
Let M2 be the space of meromorphic elliptic function on Tq with poles of order two or 
less. Then any f{z) G M2 can be decomposed in the sum of the Eisenstein functions 

/(^) = E(^2,a^2(2; - Xa\r) + Ci^aEl{z - Xajr)) + Cq, (A.34) 
a 

where 

Eci,a = 0, (A.35) 

a 

ci^a = res\^J{z), C2,a = res\^^{z - Xa)fiz), (A.36) 

and 

Co = const. part[/(z) - J2i'^2,aE2{z - Xa\r) + Ci^aEi{z - Xa\r)]. (A.37) 

a 

In particular, for a 

const. part[0(u, z — Xa)(p{—u, z — Xb)] = (p{—u, Xa — Xb)[Ei{u) — Ei{u + Xb — Xa)]. (A. 38) 
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According with (|A.34[) denote 



e2,a = E2iz - Xa), Ci^a = Ei{z - Xa) , 69 = ! (A.39) 

the basis in M2. The dual basis with respect to the integration on Tq is 

n 

f2,a = d{z- Xa)Xaiz,z), fi^a = dXa{,Z,z), fo = d{z-z){l-^Xa{,Z,z)), (A.40) 



a=l 



where Xa is the characteristic function of vicinity Ua of Xa (see ([4.25| )). 
Integrals. 

^Ci,a iEi{z - Xa\r) =^Ci^a{Xa- Xa). (A.41) 
a "^"^T a 

I E2{z - Xa\T) = -2m (A.42) 
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